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ABSTRACT 

In the framework of the internal shock scenario, we model the broadband prompt emission of 
gamma-ray bursts (GRBs) with emphasis on the GeV-TeV bands, utilizing Monte Carlo simulations 
that include various processes associated with electrons and protons accelerated to high energies. 
While inverse Compton emission from primary electrons is often dominant, different proton- induced 
mechanisms can also give rise to distinct high-energy components, such as synchrotron emission from 
protons, muons or secondary electrons/positrons injected via photomeson interactions. In some cases, 
they give rise to double spectral breaks that can serve as unique signatures of ultra-high-energy pro- 
tons. We discuss the conditions favorable for such emission, and how they are related to the production 
of ultra-high-energy cosmic rays and neutrinos in internal shocks. Ongoing and upcoming observa- 
tions by GLAST, atmospheric Cerenkov telescopes and other facilities will test these expectations and 
provide important information on the physical conditions in GRB outflows. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanisms: nonthermal 
— cosmic rays — neutrinos 



1. INTRODUCTION 

The prompt emission of gamma-ray bursts (GRBs) is 
characterized by rapid temporal variability and nonther- 
mal spectra extending to high energies, implying an ori- 
gin in ultrarelativistic outflows with bulk Lor e ntz factors 
r > 100 (see, e.g., reviews bv iPiran I [200l IMeszaros I 

2006) . In the widely discussed internal shock scenario, 
collisions among inhomogeneities within the flow lead to 
formation of shocks that convert bulk kinetic energy into 
Fermi-accelerated, power-law distributions of relativistic 
electrons, which then emit synchrotro n photons to be ob- 
serve d as the MeV-range gamma-rays (R ees fc Meszaros I 
1994). However, a number of challenges for the internal 
shock model have been pointed out concerning the radia- 
tive efficiency, low energy spectral slope, various kinds of 
luminosity correlations, etc., an d very different alterna- 
tive models have been propo sed (|Piran | [2005: Me szaros"! 
120061 : iFbx fc Meszaros II2006L and references therein). In 
order to unravel the true nature of the prompt emis- 
sion as well as to constrain important physical quanti- 
ties such as r and magnetic fields in the outflow, more 
broadband observations including the GeV-TeV bands 
are warranted. 

The physical conditions inferred for internal shocks 
indicate that protons may be Fermi-accelerated to en- 
ergies ~ 10 20 eV, making GRBs potential sources 
of the observed ultra- h igh-ene r gy co smic rays (UHE- 
CRs; iWaxman I fl995t iVietri I Il995h . To test the 
GRB origin of UHECRs and distinguish i t from 
other possibilities (jTorres fc Anchordoqui I l2004t llnoue I 

2007) , it is essential to search for characteristic, 
UHE proton-induced signatures of secondary neutral 
radiation that can be observed in coincidence with 
GRBs. Besides production of high-energy neutrinos 
(|Waxman fc Bahcall Ill997t IMeszaros fc Razzaaue 1120061 
and references therein), efficient proton acceleration may 
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induce distinctive emission components in the GeV- TeV 
bands dZhang fc Meszaros 2004: lPiran~ll2005t IMeszaros I 
120061 : iDermer fc Atovan 1120061 and references therein). 

So far, observational information on GRB GeV-TeV 
emission has been quite limited. The EGRET instru- 
ment onboard CGRO was able to detect GeV emission 
from j ust a handful of the brightest bursts dHurlev et al. I 
119941: iDingus 1 120011 iGonzalez et "al~l l2003f ) . No strong 
evidence of e mission in the TeV region has been found 
to da te fe.g. IConnaughton et al. 1 1199^: lAtkins et al. I 
120051 lAlbert et al. 1 120071: iHoran et al. I I2007T ). but this 
could be largely due to the generally high redshifts 
of GRBs and the consequent attenuation by pair pro- 
duction with extragala ctic background radiation (e.g. 
iMannheim et al. |[l996l ). 

However, significant advances are expected soon with 
the launch of GLAST 1 , with greatly improved sensitivity 
and wider field of view at GeV energies. TeV emission 
from bursts at sufficiently low redshift may eventually 
be discovered through ongoing observations with cur- 
rent Cerenkov telescopes such as H.E.S.S. 2 , VERITAS 3 , 
CANGAROO III 4 , and especially MAGIC 5 with its 50 
GeV threshold and fast slewing capabilities, as well as 
all-sky detectors such as MILAGRO 6 . 

In anticipation of the observational progress, this paper 
discusses detailed theoretical modeling of GRB prompt 
emission in the context of the internal shock scenario, 
focusing on the GeV- TeV bands. Monte Carlo tech- 
niques are employed to account for cascade processes in- 
volving photon-photon (77) pair production and Klein- 
Nishina regime Compton scattering, as well as proton- 
induced processes such as photomeson interactions and 

1 http://glast.gsfc.nasa.gov/ 

2 http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html 

3 http:/ /veritas. sao.arizona.edu/ 

4 http:/ /icrhp9. icrr.u-tokyo.ac.jp/ 

5 http://magic.mppmu.mpg.de/ 

6 http://www.lanl.gov/milagro/ 
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secondary pion, muon, electron and positron injection. 
Although various aspects of high energy emission from 
inter nal shocks have been covered in p r evious studies 
(e.g. [Papathanass iou fe Meszaros I 119961: iPilla fc Loebl 
1998HGuetta fe Granot 1120031: iPe'er fc Waxman Il2004bt 
Razzaque et al. ll2004HBaring Il2006h . few have discussed 
hadronic cascade processes in such detail. 

In fj2] our model assumptions, methods and choice of 
parameters are explained. summarizes some general 
aspects of the high-energy cutoff and inverse Compton 
emission. The effects induced by high-energy protons 
are highlighted in fj4j and the relation between GeV- 
TeV emission and UHECR and neutrino production is 
discussed in We briefly touch on the observational 
implications in Ej6j and conclude in $7] 

2. MODEL DESCRIPTION 

2.1. Model Assumptions and Numerical Methods 

In the internal shock picture, each pulse observed in 
the MeV light curve of GRBs is interpreted as emis- 
sion from shocks formed in collisio ns between mate- 
rial travelling at different velocities (|Kobavashi et al.~l 
Il997t iDaigne fc Mochkovitch 1 Il998f ). Here we do not 
deal with the dynamics of the shocks and instead con- 
centrate on the emission properties. The emitting re- 
gion for a pulse is considered to be a homogeneous 
shell expanding with T at radii R from the central en- 
gine. We adopt I = R/T for the comoving width of 
the shell, so that the p ulse timescale in th e observer 
fra me is At = R/T 2 c jlSari fc Piran I Il997t ).(see how- 
ever I Asano fc Iwamoto Il2002f ). Note that our spherically 
symmetric formulation is equally valid for a collimated 
outflow so long as the collimation angle ^> 1/T. 

Detailed modeling of the GRB spectra including the 
rapid, irregular time variability would entail consider- 
able complexity. In this work, we choose not to consider 
the time variability in earnest and assume steady state 
conditions, at least during the pulse timescale At. For 
bursts composed of multiple pulses, we also assume for 
simplicity that all pulses within a burst are similar, i.e. 
they are emitted from N shells with identical physical 
conditions. Our results are therefore to be interpreted as 
the time-averaged spectra for each burst. 

We e mploy the Monte Carlo num e rical code of 
lAsano I ([200l and |A sano fc Nagatakil (|2006D . newly 
supplemented with 77 pair production and synchrotron 
self-absorption. All photons and particles (electrons, 
positrons, protons, pions, muons) are distributed isotrop- 
ically in the shell frame and treated in the one-zone ap- 
proximation. Being mutually affected through processes 
such as photomeson interactions and inverse Compton 
(IC) scattering, the energy distributions of photons and 
particles are simulated iteratively until they converge to 
a self-consistent steady state, which is assumed to be re- 
alized within the pulse timescale. 

The energy density of accelerated electrons in the shell 
U c is a parameter that can be directly related to ob- 
servables ( §2.2|) . The magnetic field strength B is pa- 
rameterized by fs so that its energy density Ub = 
B 2 /8ir — fsUc- Electrons are injected with a power- 
law energy distribution N(j c ) oc 7 ( T Pc in the range 
7e,min < 7c < 7c,maxj where "f is the electron Lorentz 
factor in the shell frame. The minimum Lorentz factor 
7e,min is often evaluated in the literature by giving U c 



together with the total number density n e of electrons in 
the shell, w hich can be related to t he dissipated kinetic 
energy (e.g. iKobavashi et al. Ill997l ). Instead of consid- 
ering n e , here we take 7 e ,mm to be an additional parame- 
ter, the value of which can be inferred from the observed 
spectral peak energy ( §2.2p . The maximum Lorentz fac- 
tor 7c. max is where synchrotron and IC losses limit Fermi 
acceleration. However, its value is not very crucial here, 
since our choice of p c below ( ij2.2p implies that other 
factors are more important in shaping the high energy 
spectra. 

Accelerated protons with energy density U p are 
also injected with a power-law energy distribution oc 

7p Pp (7p,min < 7 P < 7 P ,max) in the shell frame. The 
maximum proton Lorentz factor 7 P!ma x is determined 
by equating t acc = j p m p c 2 /eBc, th e Fermi accelera tion 
timescale in relativistic shocks (e.g. IWaxman IIT995I ). to 
min[i cxp , ti oss ], where t cxp = R/Tc is the comoving ex- 
pansion timescale and ti oss is the energy loss timescale 
due to sy nchrotron, IC, and photomeson cooling, as de- 
scribed in |Aiano] ( 2 005) . The minimum proton Lorentz 
factor 7 P! min is expected to be of order unity in internal 
shocks with typically mildly relativistic velocities; here 
we take 7 P! min = 10, although the exact value is irrele- 
vant for the resulting spectra . 

As in lAsano fc Nagataki I (|2006l ). we utilize experi- 
mental results for the cross sections of the reactions 
P7 — > 7i7T + , pTT°, 7J7r + 7r°, and pir + Tr~ for e' < 2 GeV, 
where e' is the pho ton energy in the proton rest frame 
(Sch admand l l200l . The process pj — ► p7r°7r° is ne- 
glected due to its small cross section. For pion produc- 
tion by 717 reactions, we adopt the same cross sections 
as the respective pry channels. The inelasticity is ap- 
proximated by K = [1 — (rrip — m 2 )/s]/2, where s is the 
center-of- momentum energy squared for the pj or 717 sys- 
tem, m = TOtt and m = 2m n for single and double pion 
production, respectively, and is the pion mass. Pion 
production via pp collisions is not considered here since 
target photons always greatly outnumber protons. 

We account for the decay of pions and muons and asso- 
ciated electron/positron injection as well as synchrotron 
and IC emission from all charged particles with the meth- 
ods of I Asano I (f2005h. The full Klein -Nishina cross sec- 
tion (e.g. lBlumenthal fc Gouldlll970f ) is employed for IC 
scattering. For synchrotron radiation from very high- 
energy electrons/positrons, quantum effects can become 
important. When the classical value for the synchrotron 
photon energy e syn = j 2 heB /m c c is larger than 10% of 
the particle energy 7 e m e c 2 , w e use a pproximate emissiv- 
ity formulae following lErber I (|1966f ). The details of this 
treatment do not affect the results significantly, as such 
synchrotron photons promptly create further pairs and 
the initial information is lost in the cascade process. For 
the same reason, we also do not distinguish between the 
cascade contributions from pions and muons. 

Newly implemented here into the Monte Carlo 
code with the appropriate cross sections are 77 pair 
production and synchrotron self-absorption by elec- 
trons/positrons. The cross section for 77 pair production 
is a± — OT<?(y), where <7t is the Thomson cross section, 



(3-y i )ln\±y-2y(2-y 2 ) 
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y is given by y 2 = 1 — (2m 2 c 4 )/[£i£2(l — cos 6*)], s\ and £2 
are th e energies of the two pho tons and 9 is their incident 
angle (Bcrcstc tskii et al]|1982t ). For synchrotron absorp- 
tion of an isotropic photon field by electrons/positrons, 
the differential cross sections for true absorption and 
stimulated emission are respectively 



da s 



(7c,£o) : 
(7o,£o) : 



l h 3 



87T£q7 c M, 

c 2 h 3 , 



87T£ 



(2) 
(3) 



where £0 is the photon energy, 7 C = 7c + ^jm^c 2 , 
u c = ( 7c 2 - l)Va, U ' e = ( 7 'e - I) 1 / 2 , and P( 7c ,£ ) 
is the synchrotron power per unit photon energy 
(|Ghisellini fe Svensson |[T991f ). An accurate treatment of 
synchrotron self-absorption is necessary to determine the 
correct photon spectrum at very low energies, which in 
turn is essential for properly evaluating the photomeson 
interaction rate for UHE protons. 

We do not include pair annihilation, which can lead 
to a prominent spectral component for sufficiently high 
compactness parame ters, but only in a narro w energy 
range around Tm e c 2 (jPe'er fe Waxman~ll2004bl ). 

2.2. Constraints on Parameters 

The full set of our model parameters consists of T. 
R, N, U e , f B , 7e,min, Pc, U p and pp. For T and R, 
we consider the ranges T = 30-1000 and R = 10 13 -10 16 
cm, as generally discuss ed for internal shock models (e.g. 
iMeszaros fe Rees l2000f ). We assume a range of fs = 0.1- 
30 for the magnetic field (see below). In order to keep 
the scope of the current study tractable, some combina- 
tions of the remaining parameters are constrained so as 
to reproduce typically observed properties of the MeV, 
primary synchrotron component. 

For given values of B and T, 7 c .min is chosen so that 
the corresponding synchrotron photon energy in the ob- 
server frame £ p k = T^ 2 min he / m c c is always 300 keV (for 
GRB redshift z = 0.1, see below). The electron injection 
index is fixed to p c = 3, implying that in the fast-cooling 
conditions of internal shocks, the photon index immedi- 
ately above £ p k is /3 = — (jp + 2)/2 = —2.5, t he mean 
value measured by BATSE (jPreece et al. 1120001 ) . 

During the pulse timescale At, the fast-cooling elec- 
trons reach steady state where U c ~ Uj >e , the energy 
density of photons emitted by electrons in the shell rest 
frame. The isotropic-equivalent energy of photons from 
a single pulse is thus E s ^ e = 4:TrT 2 U^, e R 2 cAt ~ 4nR 3 U c . 
In all cases studied below, the emitted luminosity is dom- 
inated by MeV synchrotron photons, so for given R, 
U c can be related to the observable MeV pulse energy 
E s h. Hereafter U c is replaced by E B h as a parameter 
in the range 10 50 -10 52 erg. Under our assumption of 
N identical pulses constituting a burst ( q2.ip , the time- 
integrated, isotropic-equivalent photon energy for a burst 
is E to t = NE s h, which we fix to a typical value of 10 53 
erg. 

Although the proton component cannot be strongly 
constrained from existing observations, we assume U p — 
U e and p p — 2, which are necessary conditions for GRBs 
to be energetically viable as UHECR sources (|Waxman I 
1 19951 : IVietri | [l995). (However, recent observations may 



suggest larger values of U p , [J5]) The proton spectral in- 
dex p p is expected to be similar to p c at low energies 
where the particle gyroradii overlap, but this may not 
necessarily be the case at ultra-high-energies that are 
important for photomeson interactions. In particular, if 
the nonlinear back-reaction of CR pressure on the shock 
structure is significant, a concave spectral shape may re- 
sult that is much harder at high en ergies compared to low 
energies (jMalkov fc Drurvll200ll ). even though the de- 
tails are uncertain for relativistic shocks (|Baring fc Kirk I 

[Mil). 

After specifying the observables e p k, (i and E to t to 
typical values and making plausible assumptions for the 
protons, the remaining variable parameters are T, R, 
fs and _E s h- Utilizing the observable pulse timescale 
At = R/T 2 c instead of R, we choose to characterize 
our results with the set of At, E s h, T and fs- Note 
that a relation can also be made to the pulse luminosity 
/• /-,, At. 

All spectra below are shown in terms of the observed 
fluence versus photon energy, assuming a GRB redshift 
of z = 0.1. We do not include spectral attenuation 
by pair production with the extragalactic infrared back- 
grou nd, which may be justifi ed at z < 0.1 and e < 3 
TeV (jAharonian et al. M2006D . but should be more im- 
portant for higher redshifts and photon energies. The 
potent i al eff e cts of intergalactic ca s cade emission (e.g . 
Plaea 1 119951: iRazzaaue et al. 1 12004 IWang et al.l 12004 



Casanova et al. 1 120071 : iMurase et al. 1120071 1 are also ne- 
glected. 

We caution that actual GRBs are observed with con- 
siderable dispersions in £ p k, (3 and -Etot) n otwithstandin g 
a good correlation between £ p k and E tot (jAmati 1120061 ). 
Pulses within each burst can also exhibit a variety of 
properties. Such aspects need to be accounted in future, 
more comprehensive studies. 

Note that cases of Jb = Ub/U & 3> 1 can be compat- 
ible with internal shocks in a kinetic energy-dominated 
outflow, as long as the fraction of protons and electrons 
injected into the acceleration process is sufficiently small, 
and most of the outflow energy remains in the form of 
cold or thermal protons. Indeed, the typical radiative 
efficiency expected from electro ns accelerated in internal 
shock s is only a few percent ('e.glDaigne fc Mochkovitchl 
119981) (see however [Zhang et al. II2007D . so that fn as 
large as 30 may still be consistent with this picture. Even 
in magnetically-do minated flows, shocks can o ccur under 
certain conditions (Zhang fc Kobavashi | [2005). 

3. HIGH-ENERGY CUTOFF AND INVERSE COMPTON 
EMISSION 

Before proton-induced effects are addressed in detail 
in gl we discuss some generic aspects of the high- 
energy spectral cutoff that are independent of the emis- 
sion mechanism, together with the properties of GeV- 
TeV spectra in the typical case where inverse Compton 
emission from electrons dominate. 

3.1. High-Energy Cutoff 

In Fig. [TJ we show exemplary spectra for the case of 
At = 0.1 s, _E s h = 10 51 erg, and different values of T and 
fs- Above the synchrotron peak at £ p k = 300 keV, there 
are varying levels of a second high-energy component, 
here all due to inverse Compton emission. Clear spec- 
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tral cut-offs can be seen at the highest energies, above 
where pair production with low energy photons within 
the emission region strongly attenuates the spectrum. 



oc U c /Ub = fg 1 because of the Klein-Nishina effect and 
77 absorption. Note that "electrons" here include both 
primary electrons, i.e. those directly accelerated at the 
shocks, as well as additional pairs that are injected by 
77 interactions at higher energies. 



Ef(8) [erg/cm 2 ] At=0.1 s, E sl = 10 51 erg 
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Fig. 1. — Spectra for At = 0.1 s, E s ^ = 10 erg and varying 
T and f B . The thickest curve is for T = 100 and f B = 1.0. The 
medium-thick and thin curves are for T =300 and 1000, respec- 
tively, while dotted, solid and dashed each correspond to fs =0.1, 
1, and 30. 

The high-energy cutoff energy e cut should provide an 
effective probe of the bulk Lorentz factor T, as has 
been discussed previou sly (e.g. iBaring fc Harding |[l997t 
iLithwick fc Saji|[200ll). For th e case of a pure power-law 
spectrum. lAsano h Takaharal ((2003) have obtained an 

analytical expression, e cut oc T 14 / 3 i?~ h 2 ^ 3 At 4 / 3 for p c = 3 

(or e cut oc r 26 / 5 £~ 4/5 A£ 8 / 5 for p c = 5/2). Our Monte 
Carlo simulation results reveal values of e cu t that are too 
scattered to be fit well by one simple, analogous formula. 
Nevertheless, it can be approximated roughly by 

Vl00 7 Vl0 51 crg7 V Is J ' w 

with a typically strong T-dependence. Thus, together 
with the observables At and E s h (or L), measurements 
of £ cut should provide tight constraints on T. These in- 
ferences are mostly independent of the emission process 
that shape the GeV-TeV spectra, whether IC or not. 

3.2. Inverse Compton Emission 

Detecting the IC component should be crucial for prob- 
ing the magnetic field strength, especially for larger val- 
ues of T. For T =300-1000 in Fig. [TJ the strong depen- 
dence of the IC fluence on fg is apparent. For lower 
Jb and consequently higher IC fluence, the synchrotron 
fluence is somewhat suppressed due to the greater impor- 
tance of IC cooling. In contrast, for T = 100, the spectra 
for different fs =0.1-30 are almost indistinguishable, as 
e cu t occurs at too low energies for the IC component 
to be clearly discerned. Here the dominant high-energy 
component is simply the extension of the primary syn- 
chrotron emission up to e cu t ~ 0.1 GeV. 

Fig. [2] displays more details of the spectra for a 
case where IC emission makes a distinct second peak at 
e ~ 10 GeV. The IC to synchrotron peak fluence ra- 
tio is less than the simple Thomson limit expectation 
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Fig. 2. — Spectrum for a case with inverse Compton peak 
(thick solid), for At = 0.37 s, E sh = 10 51 erg, T = 300 and f B = 
0.1. Thin curves denote separately electron synchrotron (labeled 
eSY) and inverse Compton (labeled elC) components, without 77 
absorption effects. The dashed line is a power-law extrapolation of 
the MeV-range spectrum. 

These results on IC emission from internal shocks 
are broadly consist ent with previous, more approxi- 
mate studies (e.g. iPap athanassiou & Mcsz arosl 119961 : 
iPilla fe Loebl [19981 : iGuetta fc Granot 1120031 ). ~hT JAl we 
show some quantitative relations for the spectral peak 
energies and fluence ratios in our model, which might 
offer a useful consistency check of the internal shock sce- 
nario, at least within the parameter space studied here. 

4. PROTON-INDUCED HIGH-ENERGY EMISSION 

For the range of parameters covered in this study 
f q2.2p . IC emission often turns out to be the dominant 
high-energy emission mechanism (see Figs. [TJ |8j) . Never- 
theless, we find that within plausible parameter regimes, 
distinctive spectral features can emerge at GeV-TeV en- 
ergies due to characteristic processes induced by UHE 
protons. 

One potential radiative signature of UHE protons in 
GRBs i s thei r synchrotron emission, first proposed by 
IVietri I (|1997fh Fig. [3] is a case with relatively high 
magnetic fields {Jb = 30), in which proton synchrotron 
emission makes a marked contribution to the spectrum 
at e ~ 1-100 GeV. Although the spectral bump here is 
not as prominent as some examples of IC peaks ( !)3.2[) . it 
forms a clear excess above a simple extrapolation of the 
MeV-band spectrum that may be detectable by GLAST, 
MAGIC and other facilities. 

Another example with high magnetic fields (/b = 30) 
is shown in Fig. 2] Here, a notable hardening of the 
spectrum can be seen at e ^0.01-1 GeV, caused by syn- 
chrotron emission from secondary pairs injected by pho- 
tomeson interactions. Since primary electrons alone can- 
not give rise to such distinct features, this is a unique 
effect triggered by UHE protons, which has not been dis- 
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Fig. 3. — Spectrum for a case with proton synchrotron bump 
(thick solid), for At = 0.12 s, E sh = 10 50 erg, T = 300 and 
fs = 30. Thin curves denote separately electron synchrotron 
(eSY), proton synchrotron (pSY) and muon synchroton (/^SY) 
components, without 77 absorption effects. The dashed line is 
a power-law extrapolation of the MeV-range spectrum. 

cussed before for GRB prompt emission. Note that in 
order to correctly evaluate the density of low-energy tar- 
get photons for the fry process, it is imperative to include 
self-absorption effects in the electron synchrotron spec- 
trum 
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Fig. 4. — Spectrum for a case with secondary pair synchrotron 
bump (thick solid), for At = 3.3 s, E sh = 10 52 erg, T = 100 and 
Jb = 30. Thin curve denotes the synchrotron component (eSY) 
without 77 absorption effects. The dashed line is a power-law 
extrapolation of the MeV-range spectrum. 

Although such features should serve as valuable indica- 
tors of UHE protons in GRBs, it may not be easy from 
spectral measurements alone to distinguish them from 
some cases of IC emission. However, under certain con- 
ditions, more than one emission mechanism can become 
simultaneously important and lead to double spectral 
breaks, which can only occur in the presence of accel- 
erated protons. Fig. [5] is an example where the spec- 
trum hardens above a first break at ~ 0.01 GeV from 
secondary pair synchrotron emission, and then hardens 
further above a second break at ~ 0.1 GeV from IC emis- 
sion. Spectra with such double breaks may offer crucial 
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Fig. 5. — Spectrum for a case with double breaks due to 
secondary pair synchrotron and IC emission (thick solid), for 
At = 0.12 s, E Bh = 10 51 erg, T = 300 and f B = 1. Thin 
curves denote separately electron synchrotron (eSY) and IC (elC) 
components, without 77 absorption effects. The dashed lines are 
power-law extrapolations of the spectra in the ranges 1-10 MeV 
and 10-100 MeV. 

Yet a third proton-induced process that can be sig- 
nificant is synchrotron emission fr om muons injected 
by p7 interactions, first discussed bv lAsano fc Takaharal 
(2003). In Fig. O again for high magnetic fields (fs — 
30), a muon synchrotron spectral bump is eminent at 
£ ~10-100 GeV. Additionally visible in this case are sec- 
ondary pair synchrotron emission at e ~ 0.1-1 GeV, pro- 
ton synchrotron emission at e ~ 100 GeV, and even a 
minor contribution from pion synchrotron emission at 
e > 100 GeV, illustrating the spectral variety generated 
by UHE protons. While not shown here, there are other 
instances where muon synchrotron is the sole high-energy 
component (see €Bj) . 

It is important to clarify in which physical regimes of 
GRB internal shocks these proton-related emission com- 
ponents become clearly visible. Here we do not attempt 
to explore the full parameter space, but choose to map 
out certain ranges of T and fg while focusing on the fol- 
lowing two sets of At and E s ^: (1) "spiky pulse" case 
of At = 0.1 s and E sh = 10 51 erg, (similar to Figs. 
H and O, and (2) "broad pulse" case of At = 10 05 s 
and -E s h = 10 52 erg. Whenever IC, proton synchrotron 
or secondary pair synchrotron emission create distinct 
spectral features over the extrapolated MeV-range spec- 
tra, we correspondingly indicate "IC", "PS" or "SS" in 
the Jb-T plane of Figs. [7] and [HI When two compo- 
nents occur simultaneously, they are indicated together 
with a "+" sign, while black dots signify that no sepa- 
rate high-energy component is discernible. (Muon syn- 
chrotron emission does not become significant in these 
two cases, but can be evident for other At and E s h as in 
Fig. El) 

Generally speaking, we see that GeV-TeV emission 
requires sufficiently large T regardless of the emission 
mechanism to avoid 77 absorption ( §3.ip . and that larger 
/s is more conducive to proton-induced components. We 
also see that cases of multiple components can be fairly 
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and = 10 51 erg (spiky pulse case). Black dots imply no distinct 
high-energy feature. 

common. A more quantitative summary of the model 
spectra in the current study can be found in SjBj which 
may provide a guide to searches for proton-induced sig- 
nals in future observations. 

While proton-induced emission can become clearly ob- 
servable, in all the cases studied here, it does not lead 
to conspicuously separate spectral peaks as for the IC 
emission. On the other hand, such situations may be 
possible outside of the parameter restrictions we set in 
§2.21 For example, some recent observations may point 
to U p > U c if GRBs are the origin of UHECRs (©. All 
proton-related components will then be duly increased, 
as they simply scale in proportion to U p . 

Although not explicitly addressed in this work, we 
mention that time variability should also be a crucial di- 
agnostic of emission mechanisms, as each process has its 
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Fig. 8.— Same as FigO but for At = 10 5 s and E sh = 10 52 
erg (broad-pulse case). 

characteristic timescale and dependence on photon en- 
ergy. This will be an important subject for future stud- 
ies. 

5. RELATION TO ULTRA-HIGH-ENERGY COSMIC RAY 
AND NEUTRINO PRODUCTION 

We now discuss how the above results on GeV-TeV 
emission are related to the processes of UHECR and neu- 
trino production in GRB internal shocks. Since the in- 
ternal shock model entails a wide range of physical con- 
ditions by design, the circumstances most favorable for 
each process are not necessarily the same. We concen- 
trate below on some representative cases without inves- 
tigating the full model parameter space. 

The acceleration of protons to ultra-high energies is a 
necessary but not sufficient condition for GRB internal 
shocks to be significant contributors of UHECRs, since 
the particles must also escape efficiently without suffer- 
ing significant energy losses. Although a detailed de- 
scription of UHE CR esc a pe is beyond the scope of this 
paper, following lAsano I (|2005f ). we can impose a rele- 
vant constraint that the particles in question, say, pro- 
tons with energy e p > 10 19 eV, do not lose more than 
half of their energy radiatively in a comoving expan- 
sion timescale after their injection in the shell (roughly 
^loss ^ ^cxp)- Strictly speaking, this is a minimum re- 
quirement, but if it is satisfied, we may expect that the 
higher energy particles can eventually escape as the shell 
expands and both the photon density and magnetic field 
drop rapidly (note B oc i? -3 / 2 ). We can then infer 
a lower bound on the shell radius R, or equivalently 
on r for given At and E^. From our numerical re- 
sults, this criterion for efficient UHECR production is 
approximately R > 10 14 ( J B sh /10 50 erg)°- 5 (r/300)- 1 cm 
or T > 300(At/0.1s)-°- 3 (£ sh /10 51 erg) - 2 . The depen- 
dence on the magnetic field parameter fg is weak as the 
losses are mostly due to photomeson interactions rather 
than synchrotron radiation (although the latter becomes 
more important for higher energies e p > 10 20 eV). Note 
that effective particle escape via neutron conversion is 
included in this criterion and only occurs in a narrow 
parameter range near the lower limit values. 

Turning to gamma-ray emission, we saw that distinct 
GeV-TeV components mandate high values of T from 
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77 optical depth constraints, irrespective of the emis- 
sion mechanism fi i3.1|) . In fact, the above bound on 
T for UHECR production roughly matches the bound 
from gamma-rays, at least for the two exemplary cases 
of At and E s h in fj4] (regions outside the black dots in 
the /s-r plane of Figs. [7]and[8|). Therefore, the ap- 
pearance of even the IC emission may possibly indicate 
that the physical conditions are also appropriate for ef- 
ficient UHECR acceleration and escape. Of course, the 
emergence of proton-induced emission (only for high fg) 
will be most valuable as it can directly probe important 
quantities such as U p /U e and 7 Pima x- 

Detection of high-energy neutrinos is often emphasized 
as a defin itive observational test o f the GRB origin of 
UHECRs (|Halzen fc Hooper 1120021) . However, the situa- 
tion most advantangeous for neutrino production is that 
UHE protons undergo efficient photomeson interactions 
in dense radiation fields without escaping. This favors 
small values of R or V that are contrary to and almost 
m utually exclusi ve with t he UHECR criterion, as s hown 
in lAsanol (|2005f) (see also lGialis fc Pelletier Il2005f) . For 
example, the requirement that the emitted neutrino flu- 
ence > 10~ 5 erg cm -2 in the current model corresponds 
well with T < 300(Ai/0.1s)-°- 3 (£; sh /10 51 erg) - 2 , entirely 
the opposite of the UHECR bound above. (See SjC]for a 
summary of the neutrino spectra in the current model.) 
Taking this constraint at face value, we can find some 
overlap with the lowest T cases with GeV-TeV com- 
ponents in Figs. [7] and O Indeed, the pertinent pro- 
cess is found to be secondary pair synchrotron emission, 
which is generated together with neutrinos in inter- 
actions. Yet there is also a large parameter space with 
even lower T that allows copious neutrino emission but 
very little gamma-ray or UHECR production. Although 
neutrino observations will still be indispensable to verify 
that UHE proton acceleration actually occurs in GRBs, 
the bursts that emit the most neutrinos may not be the 
ones that contribute the most UHECRs. (Such remarks 
do not apply if UHECR acceleration can occur in ex- 
ternal shocks ; IVietri 1 11995): IWaxman fc BahcalTI 120001: 
Dermer Il2002f) ( see however | Gallant fe Achterberg 1119991 : 
Milosavlievic fc Nakar Il2006l regarding external forward 
shocks.) 

Thus wc find that the connection between UHECR, 
neutrino, and gamma-ray production in GRB internal 
shocks is ver y intimate, but not o ne-to-one and nontriv- 
ial (see also iDermer et al~l 12007ft . Further studies are 
warranted for a more complete understanding, but this 
point should be important to bear in mind for the re- 
spective observations. 

We remark that all of the above discussion is based 
on the assumption U p = U c ( W2.1\i . However, recent, 
post-SWIFT observations reveal the GRB redshift dis- 
tributi on to be skewed t o higher z than previously be- 
lieved (|Jakobsson et al. I I2006T ) . This may suggest that 
a larger energy budget with U p > U e may be necessary 
for GRB UHECR scenarios to remain viable, implying 
correspondingly higher gamma-ray and neutrino contri- 
butions. (Note that extr eme values su ch as U p ~ 10 3 C/ o 
have also been proposed; iTotani Ill998ft . 

6. OBSERVATIONAL IMPLICATIONS 

Here we briefly comment on the implications for exist- 
ing and future observations. 



Some EGRET-detected GRBs exhibited GeV emis- 
sion coinciding with the prompt emission, with spec- 
tra that are mostly cons istent with a n extrapola- 
tion of the MeV spectra (jDingus 1 1200 ll) . For GRB 
940217, there is some evidence of a separate high- 
energy component during the prompt phase, and per- 
haps in the delayed , hour-t imescale emiss ion as well 
(|Hurlev et al. Ill994h (see also iDermerl 12005(1 . While the 
latte r is likely to be associated with the external shock 
(e.g. iMeszaros fc Reesl 11994 iBottcher fc Dermer ||1998t 
IZhang fc Meszaros 11200 It llnoue et al. II2003I) . the former 
could possibly be related to some of the emission pro- 
cesses discussed here. More information is necessary to 
be conclusive, however. A markedly distinct component 
with a hard s pectrum above severa l MeV was seen in 
GRB 941017 ((Gonzalez et al. Il200l . but the fact that 
it varied on considerably longer timescales compared to 
the s ub-MeV emission may favor an external shock origin 
(e.g. iGranot fc Gue tta 20 01; iPe'er fc Waxman 1 |2004at 
IDermer fc Atovan I [200l iBeloborodov 1 12005ft . At any 
rate, much more detailed studies of the GeV prompt 
emission should become feasible soon after the launch 
of GLAST, which may detect some or all of the emission 
components discussed here. 

Although clear detections have yet to be achieved 
at TeV energies, the MAGIC telescope has conducted 
rapid follow-up observations for selected GRBs, in 
some cases overlapp i ng wi th the prompt emission phase 
([Albert et al. 1120061 12007ft . The obtained upper limits 
reach fluence levels of < 10 _7 erg cm~ 2 at ~ 0.1 TeV 
with integration times of several minutes, so our fiducial 
z = 0.1 burst should be readily detectable. Estimating 
the amount of intergala ctic attenuation with the base- 
line background model of lKneiske et al. I(|2004ft . MAGIC 
may be able to detect the proton synchrotron emission 
of Fig. [3] or muon synchrotron emission of Fig. [S] out to 
z < 1, and the IC emission of Fig. [2] to somewhat higher 
z, approaching the typical redshifts of GRBs. Thus the 
prospects are very promising for further observations by 
MAGIC as well as other Cerenkov telescopes such as 
H.E.S.S., VERITAS and CANGAROO III, and espe- 
cially the near-future upgraded facilities MAGIC II and 
H.E.S.S. II with their lower energy thresholds. 

Weak evidence of TeV photons coincident with GRBs 
have also been reported by some s urface detectors, e.g. 
MILAGRITO (jAtkins et al. 112000ft . However, the in- 
ferred energy fluxes are much higher than at MeV, which 
is difficult to explain in the current model framework 
unless extreme parameters are invoked, e.g. U 3> 
U e ([Totani 1 11998ft . More observations are anticipated 
for such facilities with their wide-field monitoring ca- 
pabilities, including air sh ower arrays like ARGO-YB J 
(|Di Girolamo et al. 1 12004ft and even the Pierre Auger 
Observatory (|Allard et aH 12005ft . 

7. CONCLUSIONS AND OUTLOOK 

Following the internal shock scenario and focusing on 
GeV- TeV energies, we have modelled the broadband 
spectra of GRB prompt emission through detailed Monte 
Carlo simulations including a wide variety of physical 
processes related to high-energy electrons and protons. 
Besides electron inverse Compton emission, it was shown 
that interesting proton-induced components such as pro- 
ton synchrotron, muon synchrotron and secondary pair 
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synchrotron emission can become clearly visible. Mul- 
tiple component spectra with double breaks may offer 
unique evidence of ultra-high-energy proton acceleration. 
The conditions favorable for GeV-TeV emission may also 
imply efficient UHECR acceleration and escape, but not 
necessarily strong neutrino emission. 

The observational prospects are very promising for 
GLAST, Cerenkov telescopes such as MAGIC (II), 
H.E.S.S. (II), VERITAS and CANGAROO III, as well 
as wide-field surface detector facilities. Such observa- 
tions should test the internal shock model of the prompt 
emission, provide new insights into the physics of GRB 
outflows and central engines, and probe the origin of 
UHECRs. 

Note that since we did not explicitly treat the dy- 
namics of internal shock formation, some aspects of 
our study may also be valid in more general scenarios, 
e.g. models involving magnetic energy dissipation (e.g. 
iFox fe Meszarosl I2006L and references therein), if elec- 
trons and protons can be accelerated with similar energy 
distributions. 

Our Monte Carlo simulations have also allowed de- 
tailed studies of the prompt optical emission, with inter- 
esting new results concerning both electron- and proton- 
induced spectral components. This will be the subject of 
a separate paper (Asano & Inoue, in preparation). 

More detailed and comprehensive investigations of 



the current problem should accommodate a wider 
range of parameters, including dispersions in e p k, 
and E t Qt , as well as the varie t y of pulse proper- 
ties ( e.g. I Zhang fc Meszarosl l2002t lAsano fc Kobavashl 
2003). Accounting for time variability is also an im- 
portant goal for the future. We note that such ef- 
fects may potentially smooth out some of the more 
subtler spectral features discussed here and hinder 
their observational discrimination, except for the case 
of sufficiently bright brights where time-resolved spec- 
tra can be acquired. Our detailed formalism for cal- 
culating complicated hadronic interactions and pair 
cascades should also be useful for other applica- 
tions, such as high-energy emission from the after- 
glow p hase (|Bottcher fc Dermerl 1199 8; Dcri shev et al. 



19991: ISari fc Esin . _ 
Pe'er fc Waxman 1 12005) . 



2001; Zhang & Mcsz aros I 



As this work was being comple t ed, w e became aware 
of a preprint bv iGupta fc Zhangl (|2007[ ) that addresses 
issues similar to this paper, albeit it with a simpler, an- 
alytic formulation. 



We thank F. Aharonian, Z. Bosnjak, E. Parizot, M. 
Teshima and especially F. Daigne for very valuable dis- 
cussions. 



APPENDIX 
INVERSE COMPTON SPECTRA 

When IC emission is dominant at high-energies ( §3.2|) . the IC peak energy ejc and the IC to synchrotron peak 
fluence (ef(e)) ratio may allow a useful consistency check of the internal shock model within the parameter space of 
the present study (i }2~2|) . In Fig. we summarize the dependence of these two quantities on Ai, E s h, T and Jb- 
Roughly speaking, £ic oc T 4 - 5 E~^At 15 , which is similar to the expression for e cut ( §3.1|) and sensitive to T. The peak 
fluence ratio does not vary monotonically with _E s h as it is affected by 77 absorption at large E^. The IC peak is 
suppressed for larger /b and disappears for fs = 30. 
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FlG. 9. — IC peak energy £jc (upper panels) and IC to synchrotron peak fluence ratio (lower panels), for fs = 0.1 (left) and / 
(right). Open, dotted and filled circles represent F = 100, 300 and 1000, while small, medium and large circle size correspond to 
10 51 and 10 52 erg, respectively. 
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SUMMARY OF HIGH-ENERGY SPECTRA 



A large variety of high-energy spectra are realized in the current model. Those with a clear excess above a simple 
extrapolation of the MeV-band spectrum are summarized quantitatively in the At-E s h parameter plane in Fig. 1101 
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Fig. 10. — Summary of high-energy spectral components in the At-E s ^ plane, with fs and T denoted above each panel. Dotted lines 
indicate equal luminosity at unit logarithmic intervals normalized by L51 = L/10 51 erg s — 1 . The symbols designate the relevant emission 
process; circles, rectangles, diamonds and triangles respectively signify IC, secondary pair synchrotron, muon synchrotron and proton 
synchrotron emission. Double circles are cases of IC emission with distinct second peaks. The number inside the symbol is the photon 
energy where the emission component becomes apparent in log eV units. Two overlapping symbols imply that the respective components 
occur at similar photon energies, whereas one symbol encircled by another stand for double breaks, with the corresponding break energies 
in log eV inscribed beside each symbol. 



NEUTRINO SPECTRA 

In the current model, two types of neutrino spectra can occur ( Asa no |[2005t lAsano fc Nagatakill2006| ) . When photo- 
pion cooling of protons is efficient, the neutrino number spectrum oc e~ 2 reflecting that of proton s, with a high-energy 
break at which the synchrotron cooling time of the parent pion or muon equals their lifetime (|Rachen fc Meszaros I 
Il998fh and a low energy break at which the photopion cooling time of the parent proton equals t cxp . In contrast, 
for inefficient photopion cooling, the spectrum has only one break corresponding to the latter and no portion oc s~ 2 . 
Fig. 1111 summarizes these characteristic break energies e^b when the neutrino flucncc > 10 _5 erg cm~ 2 . An analytical 
estimate for the high-energy break gives £„& oc TB -1 oc ri? 1 - 5 ^^ ' 5 oc r 4 Ai 15 i? sh ' 5 , which agrees with our numerical 
results. 
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